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Technical  part 

Summary:  For  this  project  we  have  closely  collaborated  with  Dr.  Naik  at  AFRL  to  merge  the 
computational  and  experimental  effort  to  investigate  the  biomechanics  of  Nvjp-1  protein  and  to  predict 
the  pH  and  ionic  effect  on  nanostructure  and  hydrodynamic  properties  of  the  protein  and  its 
assemblies.  The  Nvjp-1  protein  is  a  key  component  in  the  jaws  of  Nereis  virens,  a  species  of  marine 
worm.  It  contains  over  25  mol%  of  histidine,  which  is  believed  to  play  a  key  role  in  the  metal- 
coordinated  crosslinking  responsible  for  the  structural  stability  and  exceptional  mechanical 
performance  of  the  worm  jaw.  Understanding  the  nanoscale  mechanism  behind  this  crosslinking  and 
its  pathway  in  affecting  the  macroscopic  mechanical  behavior  of  the  material  is  crucial  to  develop  new 
bio-inspired  mechanomutable  materials  based  on  Nvjp-1.  Here,  we  use  a  combination  of  multiscale 
modeling  and  experimental  synthesis  to  understand  the  behavior  of  this  heterologous-expressed  protein 
from  the  nano-  to  the  macroscale.  We  have  built  a  bottom-up  molecular-based  model,  which  includes 
electronic -based  density  functional  theory  calculations,  atomistic  simulation  of  the  nanoscale 
properties  with  replica  exchange  molecular  dynamics,  and  an  elastic  network  model  for  describing  the 
macroscale  behavior  at  different  pHs.  This  multiscale  modeling  supports  the  experimental  synthesis  of 
a  photo-crosslinked  Nvjp-1  hydrogel  by  proving  both  the  nanoscale  mechanisms  and  mechanical 
behavior  predictions.  Our  theoretical  results  have  been  compared  with  experiments  and  we  find  good 
agreement  with  the  experimental  observations,  showing  that  Nvjp-1  forms  a  more  compact  structure  in 
the  presence  of  Zn2+  ions  with  a  suitable  pH  environment,  leading  to  the  formation  of  more  stable 
intra-molecular  metal-coordinated  crosslinks.  These  metal-coordinated  crosslinks  induce  nanoscale 
aggregation  of  Nvjp-1,  which  is  responsible  for  the  hydrogel  contraction  observed  in  experiments  and 
predicted  by  the  model. 


REMD  in  explicit  solvent  with 
varied  Zn  ion  concentrations 


Experimental  validation 
and  mutability 


pH  MD  in  explicit  solvent 
with  varied  pH  values: 
Full-length  and  truncated 

Nvjp-1 


Experimental  validation 
and  mutability 


Figure  1.  Illustration  of  simulation  protocol  of  ionic  and  pH  effects  on  protein  structure  and  material 
properties  study. 


Main  conclusion:  To  be  more  specific,  we  provide  a  complete  multiscale  modeling  framework  and 
experimental  synthesis  of  a  new  bio-inspired  hydrogels  based  on  Nvjp-1  protein.  The  jaws  of  Nereis 
virens  are  condensed  proteinaceous  structures  whose  mechanical  properties  are  tuned  by  coordination 
with  transition  metal  cations,  particularly  Zn2+.  Nereis  virens  is  a  marine  organism  and  the 
dehydration  during  biosynthesis  must  occur  in  an  aqueous  environment.  Since  the  removal  of  bulk 
water  from  the  newly  synthesized  protein  hydrogel  cannot  utilize  evaporation,  it  is  possible  that  the 
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contraction  behavior  observed  in  this  study  mimics  an  as  yet  unknown  natural  mechanism  for 
dehydration  of  worm  jaw  material  as  it  matures  into  a  final  functional  form.  The  contraction  of  acid- 
swollen  Nvjp-1  hydrogels  is  not  dependent  on  the  type  of  cation  present  in  the  salt.  In  fact,  the  initial 
driving  force  for  removing  bulk  water  from  the  hydrogel  is  likely  due  to  charge  screening  of 
protonated  histidine  by  the  anion.  The  type  of  cation  present  in  the  salt  seems  to  play  a  role  in  the 
degree  of  contraction  experienced  by  the  hydrogel,  where  Zn2+  treatment  consistently  yields  a  more 
compact  structure  compared  to  that  seen  with  Na+.  Our  DFT  calculations  show  a  stronger  interaction 
at  crosslinking  sites  arising  from  the  presence  of  Zn2+  ions,  for  which  a  three-coordinated  crosslinking 
is  slightly  preferred,  providing  the  most  fundamental  insight  of  the  mutable  mechanics  of  such 
biopolymer  material.  In  addition,  our  atomistic  model  of  the  Nvjp-1  protein  provides  a  better 
understanding  of  this  phenomenon,  and  predicts  that  the  Nvjp-1  has  different  mechanical  response  in 
different  solutions  with  various  ion  concentrations.  Structurally  and  mechanically  tunable  properties  of 
Nvjp-1  are  demonstrated  with  this  model  and  compared  against  the  experiments.  We  have  calculated 
SASA  of  Nvjp-1  with  the  atomistic  model,  and  have  found  that  the  SASA  of  Nvjp-1  decreases  as  the 
Zn2+/protein  ratio  increases.  We  have  also  calculated  the  radius  of  gyration  to  measure  the  size  of 
Nvjp-1,  and  the  results  show  that  a  more  compact  structure  is  formed  in  the  presence  of  Zn2+  ions. 
These  results  show  the  consistent  trend  for  solubility  as  with  the  values  measured  in  the  experiments. 
Our  coarse-grained  model  of  Nvjp-1  provides  mesoscopic  insights  of  it  mechanical  response,  which 
demonstrates  that  it  can  adapt  to  shapes  very  different  from  its  initial  states  by  generating  large 
mechanical  deformations,  providing  opportunities  to  predict  and  design  composite  materials  of  active 
mechanical  functions.  Our  results  suggest  that  metal-coordinated  crosslinks  play  a  significant  role  in 
achieving  the  characteristic  properties  of  this  protein  material. 

Our  study  shows  that  Nvjp-1  provides  an  excellent  example  of  mechanomutable  protein  material  to 
investigate.  By  integrating  multiscale  analysis  and  experimental  validation,  we  can  learn  the 
fundamental  material  science  concepts  at  different  scales  from  the  quantum  level  of  bond  formation,  to 
the  atomistic  scale  protein  folding  and  mesoscopic  scale  swelling  and  contraction  in  different  chemical 
conditions.  Our  multiscale  modeling  provides  the  most  comprehensive  understanding  of  the  formation, 
mechanical  behavior  and  the  tunable  function  of  such  protein  materials  under  changing  conditions. 
Inspired  by  Nvjp-1,  we  will  use  such  computational  modeling  tool  for  rational  polymer  material  design 
with  desired  tunable  mechanical  functions,  which  provides  the  rational  way  to  fully  program  the 
mechanical  response  of  the  material  in  different  conditions  and  is  missing  in  most  previous  studies. 

This  study  provides  fundamental  insight  into  mutable  mechanical  properties  of  Nvjp-1  and  other 
proteins  with  metal-coordinated  bonds  from  a  bottom-up  perspective.  It  also  provides  an  opportunity  to 
synthesize  biomaterials  with  mutable  properties  in  the  presence  of  metal-coordinated  crosslinks. 
Additionally,  these  pH-sensitive  materials  can  be  used  for  the  design  of  bio-inspired  adaptive  materials 
for  chemical  sensing  in  cyber-physical  systems  or  for  active  control  of  the  deformation  and  motion  of 
soft  robotics,  as  it  is  shown  in  our  large-scale  simulations  from  which  plant  coil  behavior  and  cypraea- 
like  geometries  are  generated. 
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